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Introduction 

For career astronauts and long-duration mis- 
sions, the high-energy heavy-ion component of galac- 
tic cosmic rays may cause harmful radiobiological 
effects. To assess these effects, an accurate theoret- 
ical description of the transport of high-energy nu- 
clei through spacecraft structures is being developed 
(ref. 1). An essential part of this transport theory is 
the description of the fragmentation that occurs in a 
high-energy, or relativistic, heavy-ion collision. Rel- 
ativistic heavy-ion (RHI) collisions are divided into 
peripheral, in which there is small overlap between 
colliding projectile and target nuclei, and central, in 
which there is near or complete overlap between pro- 
jectile and target. The peripheral collisions are dis- 
tinguished (ref. 2) by projectile fragments observed 
in the forward direction that are close in mass num- 
ber and velocity to that of the projectiles and by 
smaller fragments observed at large angles. In the 
much more violent central collisions, a multiplicity 
of small fragments are distributed isotropically. The 
combination of a much higher probability (ref. 3) 
for a peripheral event (approximately 90 percent) 
and the larger, higher energy peripheral fragments 
makes the description of peripheral collisions most 
important for space radiation protection. 

A previous paper (ref. 4) presented an abrasion- 
ablation T-matrix formulation for the description of 
projectile fragmentation in peripheral RHI collisions. 
The abrasion-ablation model used is illustrated by 
the schematic diagram of figure 1. Here the abrasion 
stage is represented by the interaction of the projec- 
tile P and target T, exciting these into the projec- 
tile and target prefragment states, P f and T* . To 
simplify the calculations, the removed or abraded 
nucleons from P and T in the abrasion stage are 
represented by one fireball piece R. As a further 
simplification, only the ablation stage for the projec- 
tile prefragment is considered as shown by the decay 
of P f into the final projectile fragment Z and ab- 
lated piece X. (An equivalent treatment of target 
fragmentation follows simply by interchanging roles 
of projectile and target.) The diagram of figure 1 
shows abrasion-ablation interactions to only first or- 
der. The T-matrix formalism allows for considera- 
tion of interactions to infinite order and incorporates 
ablation effects in a more fundamental way than pre- 
vious treatments (refs. 5-7). Individual abrasion and 
ablation T-matriees were shown to factor out in ref- 
erence 4, so that an optical model could be used to 
describe the abrasion stage (refs. 8-10). The main 
obstacle left in the application of this formalism is 
the calculation of the ablation T-matrix. 


In this work we use a parameterization of the 
ablation T-matrix to fit doubly differential cross sec- 
tions for projectile fragmentation in collisions of 12 C 
onto 12 C and 108 Ag at 86 MeV/nucleon, and 16 0 
onto 9 Be at 2.1 GeV/nucleon. Two parameters used 
in our calculations are the overall heights of the 
cross-sectional curves and a total width associated 
with the prefragment, or intermediate, state. Fits of 
this width allow us to make predictions of the life- 
time of the intermediate state. A width associated 
with an exponential parameterization of the ablation 
T-matrix is treated both as a parameter and as a 
temperature which is calculated as a function of the 
excitation energy of the prefragment. The fits of this 
work are compared with previous ones. In the next 
section, we review the abrasion-ablation T-matrix 
formalism; then we describe the method of evalu- 
ating Lorentz invariant and doubly differential cross 
sections; and finally we discuss our comparisons with 
experimental results and previous treatments. 

Abrasion-Ablation T-Matrix 

From reference 4 the full interaction potential 
V for the fragmentation process is separated into 
individual abrasion and ablation pieces V 1 and V 2 , 
respectively. The T-matrix expansion for abrasion- 
ablation ( AA ) is 

Tki A = Vki + VknGmV^ 

+ VknGni\ r ni 

+ V knGni V ni + V kn G ni V ni + - U) 

where G is the Green's function; subscript i indicates 
the initial state, k the final state, and n the interme- 
diate state; and an Einstein summation convention 
on triply repeated indices is implied. In reference 4, 
diagram techniques, including a time-ordering ap- 
proximation to assure causality between abrasion 
and ablation, were used to show that the abrasion- 
ablation T-matrix reduces to 

rpAA 'T'abl/nr /rtabr (o\ 

X ki ~ 1 kn Lr ni 1 ni \ Z ) 

where Tj^ and T^ v are the T-matrix expansions 
for ablation and abrasion, respectively. This is an 
extremely useful result as it allows the infinite-order 
ablation process as represented by T to be in- 
corporated into a previously developed infinite-order 
abrasion formalism. 

The T-matrix is simply related to the transition 
rate and total cross section for the fragmentation 
process. The Lorentz invariant differential cross 
section, as derived in reference 4, is most useful 
because once calculated, doubly differential, singly 



differential, and the total cross sections can be easily 
evaluated using algebraic relations and numerical 
integration. The Lorentz invariant differential cross 
section for the production of a projectile fragment Z 
is 


d 3 <j(Z) 

c 3 (dp 3 /e) z 
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where h is Planck’s constant, v is a normalization 
volume, and v the incident projectile velocity. For 
convenience we use the notation 


= 


(fajZ) 

c 3 {dp 3 /s) z 


( 4 ) 


Using the classical probability approximation 
which ignores all interference terms (discussed in 
ref. 4), 


\T^G m T t f'f « |T^‘| 2 |G m | 2 |T n t r | 2 (5) 

and expressing the abrasion T-matrix in terms of 
the total abrasion cross section a n (A), equation (3) 
becomes 

da hw (Z) = £ ^j^\T^\ 2 \G m \ 2 a n (A) (6) 

Then inserting the Green’s function through 


The subscript frame in equations (9) and (10) 
denotes that the specific quantity is evaluated in 
the desired reference frame, which is completely 
arbitrary. 


Method of Evaluation 


The calculation of Lorentz invariant differential 
cross sections for the production of a projectile frag- 
ment Z as given by equation (8) involves three main 
factors: the abrasion cross section a n (A), the abla- 
tion T-matrix Tj^\ and the Green’s function G m *. 
Calculation of each of these factors is described in 
this section. We evaluate dcr inY (Z) in the rest frame 
of the projectile as a function of the fragment’s 
4-momentum. The experimental results for the dou- 
bly differential cross sections that are used for com- 
parison are given in the lab frame. It is important 
to note (ref. 11) that the procedure for evaluation of 
d(j\ nv (Z) is to specify the fragment’s 4-momentum in 
the lab frame and to Lorentz transform this to the 
projectile frame where da\ UY (Z) is then evaluated. 

The projectile frame (fig. 2) is defined to be 
moving parallel to the lab frame with a velocity 
P = v/c, which is determined by the incident ki- 
netic energy K of the projectile per nucleon through 
(ref. 12) 


and 


1 

(l — /? 2 ) 1/2 

(11) 

K 

(12) 

= 1 + 


m N 


\G ni r = 


(s n -e t ) 2 + ( r/2)2 


( 7 ) 


where T is the total width for all intermediate states 
(projectile prefragments), we have 


where mj\f is the nucleon mass (rn^y — 939 MeV/c 2 ), 
and the speed of light c is set equal to 1 . With projec- 
tile frame quantities represented by barred symbols 
and lab frame quantities by symbols without bars, 
the needed transformation equations are 
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( 8 ) 


Equation (8) is the basis from which all numerical 
calculations are made in this report. Doubly differ- 
ential cross sections are calculated from 


d?a(Z) 

dpdU 


-1 frame 
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7 (p cos 6 — (3c) 

Abrasion Cross Sections 


(14) 

(15) 


and 


d 2 a{Z) 
de dil 


— (P c )frame^inv(^) (10) 


J frame 


Details of the calculation of abrasion cross sec- 
tions using an optical model potential approxima- 
tion to the exact nucleus-nucleus multiple scattering 
series can be found in references 8 to 10; here we 
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summarize the main results. The cross section for 
abrading n projectile nucleons is 



P{b) A P' (16) 


where b is the impact parameter vector, Ap and Apt 
are the projectile and projectile prefragment mass 
numbers, respectively, and P{ b) is the probability 
of not removing a projectile nucleon in the collision, 


P{ b) = exp [— A P a{e) 7(b)] (17) 


with 


7(b) - [2t r5(e)]- 3 / 2 | dz 3 J 
J d 3 yp P (b + z + y + Z T ) 

[l-C(y)] (18) 


x exp 


-y 
2 B(e) 


where z is the position vector of the projectile, £ 
is the collection of relative coordinates, and y is 
the two-nucleon relative position vector. Methods 
for determining the appropriate nuclear distributions 
pp and pp and constituent-averaged nucleon-nucleon 
cross sections a(e) are described in references 13 
and 14. Values for the nucleon- nucleon scattering 
slope parameter B(e) are obtained from the pa- 
rameterization, appropriate for diffractive scattering, 
given in reference 15. The Pauli correlation is 

C(y) = 0.25 exp (-4i/ 2 /lo) (19) 

where the Fermi momentum is 


K f = 1.36 fm" 1 


Results are obtained using numerical integration 
techniques. 

Since the abraded nucleons consist of protons and 
neutrons, a prescription for calculating the charge 
dispersions of the prefragments is needed in order to 
calculate final isotope cross sections. Two such pre- 
scriptions exist: a hypergeometric distribution and a 
model based on the zero-point vibrations of the gi- 
ant dipole resonance of the projectile nucleus. These 
are described in reference 16. For the calculations of 
this paper we use the hypergeometric distribution. 
Results are listed in table 1. 

Ablation T-Matrix 

The ablation T-matrix Tj^ describes the decay 
of the projectile prefragment into the final projectile 


fragment. Difficulties arise in the calculation of Tj^ 
due to lack of knowledge of intermediate state wave 
functions and of the ablation interaction. Informa- 
tion on the type of decay, such as pre-equilibrium or 
equilibrium, and accurate values for the excitation 
energy of the prefragment would be of great assis- 
tance in formulating a model for this calculation. For 
this work we looked for a simple functional parame- 
terization of T^kl which would agree well with exper- 
imental results so that some information on the abla- 
tion process could be learned. Various exponentials 
and moments of exponentials in the fragment’s en- 
ergy and momentum were studied, with consistently 
good agreement, discussed below, found only for the 
form, 


which through 


is equivalent to 




r > 2 

Pz 
2 m z 


( 20 ) 


( 21 ) 


, abl l 2 oc e Pz/ 21 ^ 2 


l^nT 


(22) 


where Kz and pz are projectile frame kinetic energy 
and momentum, respectively. The slope parameter 
and Gaussian width A are related by 


A = y/fnz~e o (23) 


The forms used for Tj^ 1 are difficult to interpret 
precisely. Gaussian shapes have been used in the 
description of projectile fragmentation in the frame- 
work of shell model, harmonic oscillator, wave func- 
tions in both a sudden approximation (ref. 17), and a 
Glauber theory abrasion- ablation model (ref. 7). A 
Gaussian width has also been related to the inter- 
nal momentum distribution of the projectile nucleus 
(ref. 18). The form given in equation (20) could also 
be interpreted as arising from an equilibrium decay 
of the prefragment. The slope parameter would then 
be the temperature of the prefragment. Because the 
physical interpretation of the form assumed in equa- 
tion (20) or (22) is not clear, we treat so both as a 
parameter and as a temperature. The temperature 
can be related to the excitation energy £* (ref. 19) of 
P f through 

e* - aTl (24) 

where a is a level density constant given by (ref. 20) 

= A Z 
8 MeV 

r F h, ^ vi f 4- V\ ^ rs ■*-\ 4- ^ i~* rv 

± lie CAeitcttiuii cucig^ ui tiic piciiA^iiieiiif 10 cia- 

sumed to be the result of two contributions: a surface 


3 


energy term s sur f due to the excess surface area of P f 
and an energy term from frictional spectator interac- 
tions (FSI) £psi which consists of abraded nucleons 
passing through P f and depositing energy. Thus 

£ — ^surf “h ^FSI 

Details of the calculation of these energies are given 
in reference 21. Note that £* must be considered 
an average energy as it depends on both the impact 
parameter of the nucleus-nucleus collision and the 
probability and trajectory of the FSI. 

Evaluation of Green’s Function 

The Green’s function G n { as given in equation (7) 
is a Lorentzian form with a full width at half maxi- 
mum T. From reference 4, F may be calculated as a 
summation over partial widths corresponding to 
transitions from all intermediate states n to a partic- 
ular final state k if the ablation T-matrix is known. 
For this paper we treat T as a parameter. Fits of T 
to experimental data allow us to predict the lifetime 
of the intermediate state through the uncertainty 
principle. 

In reference 4 the energy term in equation (7) was 
shown to be equivalent to 

e n - £i = £p> - e z - ex (27) 

which indicates a variable, intermediate, virtual res- 
onance energy Epr centered about E\ + &z- Equa- 
tion (27) must be expressed as a function of the 
4-momentum of Z in the projectile frame. In terms of 
kinetic energy, rest mass energy, and the prefragment 
excitation energy, 

£ n ~ £i = Kpi + s* + rapt 

- K z - mz - Kx ~ mx (28) 

and defining 

Q = rapt — mz — mx (29) 

we have 

e n ~ £{ = K p! + E* — K Z — Kx + Q (30) 

With conservation of mass number assumed, 

Apt = A z + A x (31) 

and with B representing binding energy, 

Q — B Z A Bx ~ B P t (32) 


For the inclusive experimental results under study, 
the exact nature of X is not known and cannot be 
specified. For example, if Ax = 2, X may consist of 
two free nucleons or a nucleus with mass number 2; 
consequently, Bx is not known. Since we consider 
light projectiles, 12 C and 16 0, where only a small 
number of nucleons are emitted in the ablation stage, 
we take 

B X = 0 (33) 

For Bpt and B z we use the approximation (ref. 19) 

- t = 8 MeV (34) 

A 

Then Q becomes 

Q = -ZA X (35) 

To obtain expressions for Kpt and Kx in terms 
of K z , we treat the ablation process, in the prefrag- 
ment frame (fig. 3), as an isotropic, two-body decay 
of P f into Z and X. In this frame, K z and Kx can 
be found in terms of rest mass energies and the ex- 
citation energy. Knowledge of K z in the projectile 
frame, where it is specified, and in the prefragment 
frame gives us the velocity of the prefragment P f rela- 
tive to the projectile P. The particle X is considered 
to move with the center-of-mass velocity of all ab- 
lated particles with a mass Axmp[. Nonrelativistic 
kinematics is used because the relative velocity be- 
tween projectile and prefragment frames is much less 
than the speed of light. (A relativistic kinematical 
treatment is given in the appendix.) 

With barred symbols representing projectile 
frame values and primed symbols prefragment frame 
values, we have by conservation of energy in the 
prefragment frame 

Ep f - E f z + E f x (36) 

which expands to 

Kpt + mpt + e* — K z + m z + Kx mx 
Using the approximation for mass terms discussed 


above (eq. (35)) and defining U by 



U = £* - 8 A x 

(37) 

we have 

U — K'z + K' x 

(38) 

or equivalently 



U = 

1 / 2 1 / 2 
-m z v z A-m x v x 

(39) 
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Applying conservation of 3-momentum in the pre- 
fragment frame 


0 - p'z + p'y 


(40) 


Similarly, 


which leads to 


VX = t>rel + v'x 


(49) 


which implies that 


Then from equation (39) 








1/2 



Kx = \ m x < 

\ 

l m z ) 

2 U 

— v z ► 

= -m z v' z 

(41) 

1 ( 

m x J 

k(^)J 


v'z ~ 


2 U 


\ mz ( 1 + sf) 


(42) 


The velocity of the prefragment as seen in the pro- 
jectile frame is equivalent to the relative velocity be- 
tween projectile and prefragment frames. Using the 
Galilean transformation equations, we have 


v Z = vz + v rel 


(43) 


or 


v'z =V 2 Z + frei + z Vrel COS <t> (44) 


where (j) is the angle between the velocity vectors v z 
and v re i. Similarly, 


Vx = vx + v rel 


(45) 


(5°) 

The expressions for Q,K P *, and Kx given by 
equations (35), (48), and (50) are used to evaluate 
the Green’s function as a function of the fragment 
variables. The kinematical model used involves many 
approximations such as assumptions on scattering 
angles of the various reaction products, neglect of 
Coulomb effects, and the assumption of simultaneous 
emission of all ablated particles. The assumption of 
one value for the excitation of P f is also an approxi- 
mation as in actuality a particular prefragment may 
have various values for its excitation energy. The 
distribution of values for s* is expected to have a 
maximum of about 100 MeV for the prefragments 
studied, and thus any error in choosing s* introduces 
only a small error in equations (48) and (50). These 
approximations are necessary in order to compare 
with experimental results, but eventually a quantum 
mechanical calculation of the energy terms must be 
made. 


v'x 2 = i>x + Vrel + 2 V* Vrel COS <* (46) 

where a is the angle between the velocity vectors 
\X and v re j. There is insufficient information within 
our model to determine the angles (j) and a except 
for the special case when the fragment Z is detected 
at a scattering angle 6 of 0°. Here we may assume 
that P f and Z are moving parallel to each other in 
both lab and projectile frames, so that <fi = 0° and 
a — 180°. For the experimental results for 16 0 on 
9 Be at 2.1 GeV/nucleon, all fragments were detected 
at $iab — 0°. The experimental results for 12 C on 
12 C and 12 C on nat Ag are for fragments detected at 
$l a b = 4° and 0] a b = 3°, respectively. In all cases 
we assume that 0 = 0° and a = 180°. With this 
approximation we find that 


u re i = v P f = v z -v z 


(47) 


and thus 


K P , = ^ m pt 


1 2 


u z - 


2 U 




+ m z \ 
™<X J 


(48) 


Results and Discussion 

As a first application of the abrasion-ablation 
T-matrix formalism, we have made two- and three- 
parameter fits to experimental results of d 2 a/dpdQ 
for projectile fragments in 16 0 on 9 Be collisions at 2.1 
GeV/nucleon (ref. 22) and of d?a/dedVt for projectile 
fragments in 12 C on 12 C (ref. 23) and 12 C on nat Ag 
(ref. 24) collisions at 86 MeV/nucleon. For the 12 C 
on nat Ag reaction, we assume a target of 108 Ag. The 
normalization of the experimental curves is reported 
only for 12 C on 12 C with arbitrary units used in 
the other results. All the experimental curves have 
simple bell shapes and were fit to good agreement 
by our two-parameter model as shown in figures 4 
to 6. A nearly identical fit is made with the three- 
parameter model discussed below. The excitation 
energies and corresponding temperatures used are 
listed in table 2. In our kinematical model, there 
was insufficient excitation energy for a prefragment 
state with A — 15 for 10 C and with A = 11 for 7 Be. 
Table 3 lists the fitted values of the total width T 
for each projectile fragment versus T n of the two- 
parameter model, and table 4 lists fitted values of T 
versus so the slope parameter of our three-parameter 
model. The values of £q and the corresponding values 
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of r listed all produce a fit nearly identical to our two- 
parameter fit shown in figures 4 to 6. We were unable 
to obtain a three-parameter fit at £o = 2 for n C in 
the 16 0 on 9 Be reaction and for 10 B and 7 Be in the 
12 C on 12 C reaction. 

The fits were made to reproduce the experi- 
mental curves in their central regions. Reasonable 
agreement is found in the high- and low-momentum 
tails for the C fragments of 16 0, the low-momentum 
tail of 12 C (fig. 4(c)), and high-momentum tail of 10 C 
(fig. 4(e)) being exceptions. Some deviation is seen 
in general for the high- and low-energy tails of the 
lower energy reactions (figs. 5 and 6). The peaks of 
all curves occur at a velocity slightly less than the 
projectile velocity. This phenomenon has been ex- 
plained previously as the result of frictional (ref. 25) 
and binding energy (ref. 26) effects. 

The fitted values for the total width T are the 
same order of magnitude for all fragments stud- 
ied and through the uncertainty principle indicate 
an intermediate state lifetime r on the order of 
10“ 19 to 10“ 2 ° s. The widths for 10 B and 7 Be frag- 
ments from 12 C targets and from nat Ag targets differ 
by a factor of about 4, with a larger width and hence 
shorter lifetime for the smaller target 12 C. This com- 
parison is made at slightly different scattering angles, 
4° for 12 C and 3° for nat Ag, so that it is difficult to 
conclude that any target dependence has been found. 
The similarity in values of T at the two energies stud- 
ied can be attributed to the comparable size of the 
projectiles and the similarity in excitation energies of 
all prefragments. 

In considering the values for r predicted from our 
fits to T, we consider two other characteristic times: 
r\ ss 10“ 23 s (the interaction time for the abrasion 
stage) and r 2 ~ 10“ 21 s (the typical period for 
nucleon motion in the nucleus). A rough indication 
of an equilibrium decay or compound nucleus decay 
of an excited nucleus is that the decay time be much 
longer than T 2 . Our predicted values for r satisfy 
this criterion and thus give some justification to the 
use of evaporation codes for the treatment of the 
ablation stage, although the mode of formation of 
a compound nucleus state is very different from the 
violent nucleon knockout of abrasion. The values of 
r also indicate that r T\ which is consistent with 
the use of the time-ordering approximation (ref. 4) in 
the formulation of the abrasion-ablation T-matrix. 

From table 4, we see that as sq is increased, T 
decreases, indicating a longer intermediate state life- 
time for larger so. This seems to contradict the inter- 
pretation of sq as a temperature for which we would 
expect a shorter lifetime for higher temperature. As 


so is increased past 10 MeV, the fitted values for 
T converge to constant values. Here the exponen- 
tial parameterization of Tj^ becomes very flat and 
the Green’s function and kinematical factors domi- 
nate the cross-sectional distributions. In the central 
regions of the distributions, agreement is good for 
these higher values of sq, but agreement in the tail 
regions decreases. 

The experimental data with which we compared 
our model have been parameterized previously using 
a simple Gaussian form (refs. 22 to 24) 

_«2/2A 2 

e F 

with the width A being the focus of much investiga- 
tion (refs. 7, 17, 27 to 31). Values for A have been 
derived in one-step models of quick fragmentation 
(refs. 17 and 18) and of excitation followed by statis- 
tical decay (refs. 18 and 27). In reference 7, values 
for A were predicted in an abrasion-ablation model 
with ablation representing a correction to values for 
A derived from the abrasion step. 

In figure 7 we show fits to the data for n C 
from the 16 0 on 9 Be reaction and for 7 Be from 
the 12 C on 12 C reaction. Shown are fits from our 
two-parameter model, from the Gaussian form, and 
from the abrasion-ablation T-matrix formalism with 
|T^ l | 2 = 1. All three fits reproduce the central re- 
gions and widths of the curves, with the first two 
being nearly identical and the fit with constant Tj^ 
deviating in the tail regions. These fits indicate that 
the widths of the experimental distributions do not 
give a strong enough indication to establish the cor- 
rect model with which to describe the fragmentation 
process. If we denote the actual width of the distri- 
bution by re, then in the Gaussian models, w equals 
A; and in the abrasion-ablation T-matrix formal- 
ism, w is a function of T n and I\ The temperature 
T n and Gaussian width A may be related through 
equation (23). In reference 18, the relationship be- 
tween A and a temperature is also discussed. In the 
abrasion-ablation T-matrix formalism, deviation of 
predictions of temperature or A from experimental 
values for w in the one-step models is explained by 
the inclusion of the Lorentzian form that arises nat- 
urally through the Green’s function that propagates 
between the abrasion and ablation stages. 

Concluding Remarks 

An exponential parameterization of the ablation 
T-matrix leads to good agreement with experimental 
results for doubly differential cross sections in frag- 
mentation reactions involving light projectiles. Fits 
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to the total width for the ablation stage indicate a 
quick decay process, but an equilibrium decay. The 
importance of prefragment excitation energies is al- 
ready apparent at the level of calculation of this 
work. Explicit calculation of the ablation T-matrix 


must now be conducted in order to make further 
progress. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
December 8, 1986 
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Appendix 


Similarly we find 


Relativistic Kinematical Model 

In this appendix we give the relativistic kinemat- 
ical model for the calculation of energy terms in the 
Green’s function of equation (7), which is analogous 
to the nonrelativistic model presented in the text. 
Relativistic effects increase in importance for higher 
values of excitation energy. For the two-body decay 
of P ( into Z and X, assumed in the prefragment rest 
frame, we have the following 4-momentum 


Pp/ — (£p/, 0 ) (Al) 

Pz = {£z,Pz) (A2) 

Px = {£x, Px) (A3) 

Through conservation we have 

P Z = Pp> - Px (A4) 




(mpi + £*) 2 + m 2 z — m 2 x 
2 (mp> + £*) 


(A9) 


As discussed above, we take the prefragment’s 
velocity, as observed in the projectile frame, to be the 
velocity of the prefragment relative to the projectile, 
which is found through determination ot s z in both 
frames. By defining the lab and projectile frames 
as parallel, we are not free to choose the relative 
orientation of prefragment frame with respect to 
projectile frame (fig. 3). We therefore take, as the 
Lorentz transformation matrix between these frames, 
a boost in an arbitrary direction. With A(/3) the 
transformation matrix, we have 

P = A{(3) P (A10) 


Taking the scalar product of Pz with itself we find 
m z — e 2 p , + ra\ — 2s ptSx ( A5) 


where P and P are projectile and prefragment frame 
4-momentum, respectively. The matrix A(/3) is given 
by equation 11.98 of Jackson (ref. 32). Then from 


Then solving for ex , 




s 2 p , + rn 2 x — m 2 z 
2s pt 


Pp> = ^rel = Pz ~ Pz (All) 


< A6 > and 


and inserting 
we have 


Spt = nipt + s* 


£X 


(mp/ + £*) 2 — m\r — m| 

2 (mpt + £*) 


(A7) 

(A8) 


Px = Px - Pre\ (A12) 

we can formally solve for the energies Spr and Ex- 
Difficulties arise in that values for the parame- 
ters of the transformation matrix A(/3) cannot be 
determined. 
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Symbols 

a level density constant, MeV -1 

A mass number 

B binding energy, MeV 

c speed of light 

C'(y) Pauli correlation 

G n { Green’s function matrix 

/(b) defined in equation (18) 

K kinetic energy, MeV 

m mass 

p 3-momentum 

p momentum 

P projectile (see fig. 1) 

P f projectile prefragment (see fig. 1) 

R fireball piece (see fig. 1) 

T target (see fig. 1) 

T' target prefragment (see fig. 1) 

^kn ablation T-matrix 

T j^ r abrasion T-matrix 

rAA abrasion-ablation T-matrix 

T n temperature, MeV 

v velocity vector 

v re l relative velocity vector between projectile 

and prefragment frames 


V 1 abrasion potential 

V 2 ablation potential 

X ablated piece (see fig. 1) 

Z final projectile fragment (see fig. 1) 

a angle between velocity vectors \x and 

Vrel 

/ 3 relative velocity vector between lab frame 

and projectile frame 

T total width for all intermediate states 

A Gaussian width 

£ energy, MeV 

£o slope parameter 

e* excitation energy, MeV 

9 scattering angle 

0 cross section 

a n (A) total abrasion cross section, mb 

(j) angle between velocity vectors v z and 

Vrel (see eq. (44)) 

H solid angle, sr 

Subscripts: 

1 initial state 

k final state 

n intermediate state 

A bar over a symbol indicates the projectile frame. 
A prime (other than P f ) indicates the prefragment 
frame. See figures 2 and 3. 
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Table 2. Excitation Energies and Temperatures 


(a) 16 0 + 9 Be at 2.1 GeV/nucleon 


Apt 

£* , MeV 

T n , Mev 

15 

35.6 

16 . 8 /v^i 

14 

39.2 

17.7/ s/j^z 

13 

42.7 

18.5/v^ 

12 

46.1 

19.2/v^i 

11 

47.3 

19.4/v^i 


(b) 12 C + 12 C at 86 MeV/nucleon 


A pt 

e*, MeV 

T n , MeV 

11 

30 

15.5 / \fKz 

10 

38 

17.4/ 

9 

46 

19.2/ \J~Kz 

8 

49 

19.8 lyfXi 


(c) 12 C + 108 Ag at 86 MeV/nucleon 


A pt 

e*, MeV 

T n , MeV 

11 

25.0 

14.1/v^i 

10 

35.0 

16.7/ 

9 

41.4 

18.2/vO^ 

8 

42.6 

18.5/ 





Table 3. Fitted Values for Total Width T Using Calculated Value of T n in Two-Parameter Model 


(a) 16 0 -F 9 Be at 2.1 GeV/nucleon with 9^ = 0° 


T, MeV for — 


I3 C 


n c 

io c 

0.018 

0.014 

0.03 

0.048 

0.027 


(b) 12 C + 12 C at 86 MeV/nucleon with 6^ = 4° 


T, MeV for — 

1U B 

7 Be 

0.12 

0.09 


(c) 12 C + 109 Ag at 86 MeV/nucleon with 0i a b = 3° 


T, MeV for — 

TOg 

7 Be 

0.03 

0.02 
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Table 4. Fitted Values for Total Width T as a Function of Slope Parameter sq in Three-Parameter Model 


(a) 16 0 + 9 Be at 2.1 GeV/nucleon with 0 lab = 0 ° 



T, MeV for — 

So, MeV 

14 c 

13 C 

12 C 

U C 

10 c 

2 

0.04 

0.03 

0.2 

(a) 

0.04 

3 

.025 

.02 

.05 

0.1 

.03 

5 

.02 

.015 

.03 

.04 

.02 

7 

.015 

.012 

.025 

.03 

.015 

10 

.015 

.011 

.02 

.03 

.015 


(b) 12 C + 12 C at 86 MeV/nucleon with 0 iab = 4° 


£ 0 i MeV 

r, MeV for — 

10 B 

7 Be 

2 

(a) 

(a) 

3 

0.5 

1.0 

5 

.1 

.15 

7 

.08 

.08 

10 

.07 

.06 


(c) 12 C + 108 Ag at 86 MeV/nucleon with 0i ab — 3° 


£o, MeV 

T, MeV for — 

10 B 

7 Be 

2 

o.fo 

0.07 

3 

.05 

.04 

5 

.03 

.03 

7 

.03 

.02 

10 

.02 

.02 


a Could not be fit with listed value of Sq. 
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Figure 2. Relative orientation of lab and projectile frames. 



Figure 3. Orientation of prefragment frame. 




Momentum, p, GeV/c 
(a) For 14 C production. 

Figure 4. Doubly differential cross section for fragment production at 0° in reaction of 16 0 on 9 B 
at 2.1 GeV/nucleon. 


■(j/dpdQ, arbitrary units d 2 o / dp dVl. arbitrary units 



Momentum, p, GeV/c 

(b) For 13 C production. 



Momentum, p, GeV/c 

(c) For 12 C production. 
Figure 4. Continued. 
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Experiment (ref. 22) 
Two-parameter fit 
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